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We discuss stationary aspects of a set of driven lattice gases in which hard-core particles with 
spatial extent, covering more than one lattice site, diffuse and reconstruct in one dimension under 
nearest-neighbor interactions. As in the uncoupled case [M. Barma et al, J. Phys. Condens. Matter 
19, 065112 (2007) ], the dynamics of the phase space breaks up into an exponentially large number of 
mutually disconnected sectors labeled by a non-local construct, the irreducible string. Depending on 
whether the particle couplings are taken attractive or repulsive, simulations in most of the studied 
sectors show that both steady state currents and pair correlations behave quite differently at low 
temperature regimes. For repulsive interactions an order-by-disorder transition is suggested. 
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I. INTRODUCTION 



Since a general theoretical framework for studying nonequilibrium phenomena yet remains elusive, our understand- 
ing of the subject partly has to resort to studies of specific and seemingly simple stochastic models. One of the most 
investigated ones in that context is a driven lattice gas (DLG) involving hard-core particle diffusion under an external 
bulk field which biases the particle flow along one of the lattice axes. Introduced by Katz, Lebowitz and Spohn [l| 
in part with the aim of studying the physics of fast ionic conductors it has triggered a great deal of research 
for almost three decades [1, 0] ■ In the high field limit this system describes an asymmetric simple exclusion process 
(ASEP) which aheady in one-dimension (ID) under suitable boundary conditions has encountered applications 

as diverse as protein synthesis [t^ , inhomogeneous interface growth ^ , and vehicular traffic ^ . Notwithstanding the 
deceptive simplicity of most DLG versions, actually slight modifications of kinetic Ising models [10], the constantly 
maintained bias results in a net dissipative current (if permitted by boundary conditions), so the emerging steady 
state (SS) distributions are nonequilibrium ones. 

As part of the ongoing effort in this context, here we investigate numerically stationary aspects of DLG with 
extended objects which in turn can dissociate and reconstruct themselves in a ID periodic lattice. Although less 
frequently studied, exclusion processes with spatially extended particles date back to the work of MacDonald, Gibbs, 
and Pipkin [TT| . who introduced this concept as a simple setting to understand the dynamics of protein synthesis. In 
that terminology, each lattice site denotes a codon on the messenger RNA, while large particles stand for ribosomes 
which, covering several codons, move through them stepwise and thereby produce the proteiii]lll[. Subsequently, this 
and other issues related to diffusion of extended objects were revisited in various studies 0, [l2| - [l8| . 

In common with some of these latter, the processes here considered involve hard-core composite particles, hereafter 
termed fc-mers, which occupy k consecutive locations and diffuse by one lattice site in the presence of both an external 
drive and other fragments of length I < k. Besides, we include nearest-neighbor (NN) interactions between individual 
particles (hard-core monomers) , and allow for fc-mer dissociation (l3l - [l5| in the course of their casual encounters with 
the otherwise non-diffusing fragments, e.g. ••o»o ^ o»»»o o«o»», say for dimers approaching a monomer. 
Thus, although the dynamics preserves the total number of A:-mers, their indentities as well as those of the fragments 
collided are not retained, these being instead recomposed throughout the process. At finite drifts and couplings the 
system evolves to a nontrivial SS measure characterized by a macroscopic current, a central quantity of interest which 
in the terminology of protein synthesis corresponds to the stationary protein production rate 0. When monomers 
are uncoupled (beyond their excluded volumes), all SS configurations are equally likely under periodic boundary 
conditions (PBC) but the full dynamics and stationary correlations are quite involved [isl. [Tsj. 

Interestingly, and irrespective of the monomer couplings, ergodicity is strongly broken as a result of the presence 
of the aforementioned fragments. As these do not diffuse explicitly (but only through dissociation with fc-mers), they 
break up the phase space into mutually disjoint and dynamical invariant subspaces or 'sectors' whose number turns 
out to grow exponentially with the lattice size [isl-fisl Il9|. Thus, the SS current and correlations are not unique but 
rather vary from one sector to another, ultimately depending on the initial distribution of fragments. In that regard, 
to characterize the full partitioning of the phase space here we follow the ideas given in Refs.[l3l-[l5| for the case 
of dimers, and introduce a set of fictitious extended particles whose order along a non-local construct -namely the 
'irreducible string' [19]- comes out to be the actual invariant of the motion. This also enables us to obtain saturation 
currents of generic sectors by means of a correspondence to ASEP systems [T^ . When normalized to those saturation 
values, our simulations indicate that as long as monomer interactions are held attractive the currents of most studied 
sectors can be made to collapse into a single universal curve. Moreover, this latter can also be fitted in terms of 
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mean-field DLG currents [J, |20[ upon using the ASEP densities associated to each sector. By contrast, for repulsive 
interactions such normalized currents turn out to be sector dependent and no universality can be constructed. On 
a mesoscopic level of description, also different features show up depending on whether the particle couplings are 
taken attractive or repulsive. Although in either case most sectors bear highly degenerate ground states, based on the 
behavior of both the structure factor and correlation length at low temperature regimes, we suggest that for repulsive 
couplings thermal fluctuations appear to lift part of this degeneracy and cause an order-by-disorder transition pll |. 

The layout of this work is organized as follows. In Sec. II we define the basic kinetic steps and transition probability 
rates of these reconstituting processes. We then recast the dynamics in terms of new extended particles which readily 
evidence the appearance of an exponentially growing number of disconnected subspaces. Also, for large drives these 
new particles are helpful to characterize the mentioned analogy between reconstructing DLG and ASEP systems. 
Guided by these developments, simulations for dimers and trimers are discussed in Sec. Ill where we examine SS 
currents and pair correlations in several sectors under various situations. We close with Sec. IV which contains a 
recapitulation along with brief remarks on open issues and possible extensions of this work. 



II. DIFFUSION OF COMPOSITE PARTICLES 



The microscopic particle model we consider is a ring of L sites each of which may be singly occupied (occupation 
number n — 1), or empty {n — Q). The particles behave as if they were positive ions in relation to a uniform electric 
field £, while in turn are coupled effectively either by NN attractive interactions (J > 0), or NN repulsive ones (J < 0), 
via an Ising Hamiltonian H — — 4 niHii^i. Let us first describe the basic kinetic steps which take place and 
then carry on with the definition of their corresponding rates. The system evolves stochastically under a particle 
conserving dynamics involving just fc-mer shifts in single lattice units, i.e. 

1_^0 ^ 1_^, (1) 

k k 

the motion being biased in the direction of the field. Here, monomers and groups or fragments of j-adjacent particles 
with j < k can not diffuse explicitly but since the identity of fc-mers is impermanent, they are ultimately allowed to 
in a series of steps. For instance, in the sequence 

1 ... 1 ^ L .1 1 ... 1 ^ 1^ 1 ... 1 , (2) 

k j j k j k 

the rightmost group of j particles can hop fc-sites to the left and vice versa. The key issue to emphasize is that both 
fc-mers and fragments can dissociate and reconstitute without restrains throughout the process, so they do not retain 
their indentity (except in particular situations, as we shall see below). 

As for transition rates between two particle configurations C and C, we take up the common Kawasaki transitions 
(I}[I3{AH + u£ k)] with <j){x) = 2(l-He^)-\ and AH = H{C')-H{C) Here, the term uffc denotes the work 

done by the bulk field £ during a fc-mer jump (u = ±1), whereas /3 = 1/T stands for the usual inverse temperature 
(henceforth the Boltzmann constant is set equal to 1). The different situations encompassed by these rules along with 
the rates associated to them are easily visualized in Table U (where thereafter Ek = (5 S k). In the absence of drive 
the SS distribution is of course proportional to e.~^'y'^^l^ but owing to PBC, when £" 7^ these rates no longer satisfy 
detailed balance (23| and the exact form of the nonequilibrium SS distribution is unknown. 



Process Rate (^) 



i [1 ±tanh(£'fc/2)] 
i [l±tanh(Sfc/2)] 
i[l±tanh(£fc/2-2/3J)] 
i [l±tanh(i;fc/2-h2/3J)] 



TABLE L Transition rates of driven fc-mer processes in Id. The symbols • , o, and | | denote respectively occupied, vacant, 
and fc-mer locations, whereas upper and lower signs stand in turn for forwards and backwards •<— hopping rates. 
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This set of driven and reconstituting lattice gases (DRLG) may be viewed as an extended and interacting version 
of a dimer model recently studied in Ref. 15| . Alike its non-driven (equilibrium) predecessors [H, , the dynamic 



here considered splits up the phase space of configurations into an unusual large number of invariant sectors, actually 
growing exponentially with the system size (see Sec. II B). Apart from particle conservation within each of the k 
sublattices (say for integer L/k), the full partitionin g of the phase space can be understood in terms of a nonlocal 
construct known as the irreducible string (IS) [T54l5n i9l|. This latter is an invariant of the stochastic motion and in 
turn provides a convenient label for each sector. To explain further this idea we recur to an equivalent representation 
of these processes using a set of composite characters or new 'particles' Aq ,Ai , . . . , A^-i , A^ constructed as 

= 0, 

Aj = l^.AQ , l<j<k, (3) 




Ak ^ 

The Ak or fc-mer movements and their recompositions can then be thought of as character exchanges of the form 

Ak Aj ^ AjAk, 0<j<k, (4) 

the fc-mer identity here being preserved only by Aq, whereas exchanges not involving Ak remain disabled, i.e. Ai Aj 
do not swap their positions if «,j ^ k. For example, in this representation the steps referred to in ([2]) now become 
Ak Aq Aj ^ Aq Ak Aj ^ Aq Aj Ak ■ This bears some resemblance to a driven process of several particles species 
introduced in the context of ID phase separation [2^. However, in those systems all NN species are exchangeable 
[25} whereas in this mapping note that the ordering of characters Aq, Ai, ■ ■ ■ ,Ak-i set by the initial conditions is 
conserved throughout all subsequent times, modulo eventual interposition of one or more adjacent fc-mers between ^'s. 
Thus, the invariant IS of a given sector simply refers to the sequence emerged after deletion of all fc-mers or 'reducible' 
characters appearing in any configuration of that sector. This operation results in a unique sequence irrespective of 
the order of deletion, so the correspondence between original monomer and character configurations is one-to-one. 

For this new representation we can therefore think of an equivalent Haniiltonian Ji of hard-core particle species Ai 
of fixed concentration Ni/N, defined on a ring of iV = J^i sites with NN interactions. In terms of the occupation 

(i) 

numbers n\ of these new particle classes, the related Hamiltonian can be written down as 

j#0 A X 

up to a sector dependent constant (here, if n^^^ = 1 at a given location A, then n^-'-* = 0, Vj ^ i). Hence, the Kawasaki 
rates of the biased exchanges referred to in ^ now depend on the class of surrounding particles, namely 

0[/3(A«-£fc)] 

Ai Ak Aj Ai ^ A, Aj Ak Ai , j^k, (6) 

0[-^(A«-£fe)] 

with 4>{x) defined as above, and energy changes given by 

AH ^ / if ^ 7^ or ? = j = , i 7^ fc , 
I 4 J otherwise . 

In particular, it follows that all IS sectors having Nq = conserve the internal energy throughout. 

It is worth pointing out that for fc = 1 the form of Eqs. ([S][7]) just describes the usual monomer DLG. Also, notice 
that the dynamics of this latter is formally analogous to that of the non reconstructing or null sector [Aq]^, that is, 
an IS of length £ — Nq. In that case Eq. ([5]) reduces to the usual Ising Hamiltonian, and the processes of © just 
involve Ak-Ao ('particle- hole') exchanges. In passing, we mention that it is actually the non-interacting version of 
this sector the one which was studied in connection to the protein dynamics referred to in Sec. I, and the one whose 
space-time correlations were recently investigated in Ref. |18| . 



A. The ASEP limit 



For iVo = as well as in the large field or saturation regime \£\ ^ | J|/fc of generic sectors, clearly the above processes 
are isomorphic to an ASEP system in which fc-mers play the role of non-interacting (but hard-core) particles hopping 
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through indistinguishable Ai vacancies {i ^ k). More specificaUy, given an IS of length C = X^i^^fcC* + 1) -^i measured 
in the original DRLG spacings, then this regime amounts to a problem of Nk = {L — C)/k ASEP particles hopping 
with biased probabilities ^(1 ± tanh ■^) through N = Ni lattice sites with particle density 

From these correspondences we can readily obtain the sublattice currents of the original DRLG in their saturation 
regimes. Since in our ASEP limit all SS configurations are equally likely (because of PBC), evidently the probability 
of finding an Ak particle followed by an Ai vacancy or vice versa, is just p^^s^p (1 ~ Pasep) j ^ turn proportional to 
the ASEP current. In the DRLG representation this is related to the probability of finding a /c-mer 'head' ('tail') 
- i.e. a rightmost (leftmost) fc-mcr unit - on a given sublattice site, times the probability of finding a in the next 
(previous) sublattice location. However, the former event occurs with probability Nk/{L/k) = Pj^g^p k N/L , whereas 
that of the latter, cx (1 — P^sep)^ must be normalized by the fraction of O's of the sublattice in question. So if pi 
denotes the monomer density of sublattice A^, such fraction is therefore calculated as 

hence, for \£\ ^ \ J\/k the saturation current of A^ finally reduces to 

J1:1 = Pasep (1 " P^) tanh {Ek/2) , (10) 

while also holding V£ so long as Nq — 0. In general, these saturation currents depend on the particular distribution 
of characters in the IS. However for periodic sectors, i.e. strings formed by repeating a unit sequence of A' s, these 
limiting values are just rational functions of sublattice densities (cf. Table HIl below). In particular and for ulterior 
comparisons, in the null sector [Aq]^ referred to above all sublattices share a common saturation current 

_ p(l-p) tanh(ii;fc/2) 

and a particle density p = 1 — C/L. 

Although in Sec. HI we shall restrict ourselves to time independent SS aspects, such as currents and pair correlations, 
let us briefly mention here that density fluctuations in the stationary ASEP move through the system as kinematic 
waves that ultimately take over the asymptotic behavior at large times. In our DRLG model this corresponds 

to k sublattice wave velocities Vi = djj ^] /dpj which, for periodic strings, can vanish at a common critical length 
Cc- Hence, as suggested in Refs. [isl. Il8j| when approaching such conditions there may well be a crossover from an 
exponential relaxation of density fluctuations to a slow Kardar-Parisi-Zhang dynamics 6] for which the former would 
decay as t~^^^. In particular, from Eq. pTt it follows that in the null sector this should occur for pc Vk/ (1 + Vk) 



B. Growth of invariant sectors 



Before proceeding to the simulation of these processes at finite temperatures and fields, we pause to digress about 
the exponential growth of disjoint sectors with the length of their IS's. Here we follow and extend slightly the 
recursive procedure discussed in Ref. [l^ for the case of dimers. For simplicity, and solely for the purpose of avoiding 
the overcount of strings related each other by a cyclic permutation of characters, we will assume open boundary 
conditions throughout this subsection. 

Let Fc (1) and Fc (0) denote the number of IS's of length C whose first bit is 1 and respectively. Thus, the total 
number of C- sectors we want to evaluate is just A^^ — Fc (0) + Fc (1). As there are no A^ characters in these strings, 
then these quantities must satisfy the recursion relations 

^£(0) = Fc-i{0)+Fc-iil), 

Fc{l) = Fc-iiO)+Fc-iiAi) + ... + Fc-iiAk-2): (12) 

with Fc (Ai) denoting in turn the number of IS's of length C whose first character is^^, l<i<fc— 1 (note that 
such construct would not be well-defined for PBC). On the other hand, by definition these latter quantities should 
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also be related recursively as 

FciAi) = Fc-iiO), 

■■ (13) 
Fc (Ak-i) = Fc-i {Ak-2) = ■■■ = Fc-(k-i) (0) • 

As a byproduct of these relations, it follows that Nc = 2Fc (0) — Fc~k (0), thus it is sufficient to focus attention on 
Fc (0). Inserting the above equations in Eq. (|12l) we readily obtain a linear recursion for this quantity, namely the 
Fibonacci-like relation 

fc 

Fc{0)=Y,Fc-^{0) , Ok, (14) 

1=1 

whose general solution is bound up to the zi , . . . ,Zk zeros of the polynomial [l^l 

fc 

Pk{x)=J2^'~'- ■ (15) 

i=l 

Thus, Eq. reduces to the exponential form Fc (0) = J2^=i ^'^ii^ ^-coefficients that are in turn evaluated by 
fitting linearly the k boundary terms i^i(O), ... ,Ffc(0), e.g. Fi(0) = F2(0) = 1 for dimers, and Fi(0) = ^2(0) = 
^^3(0) = 1 in the case of trimers. Specifically, for these two latter situations, which we shall discuss in detail for 
some sectors later on, it turns out that for large C (say comparable to L), the total number of £-strings grows as 

f [i(l + %/5)]^~1.618^, for fc = 2, 
Nco,\ (16) 

[ [i {1 + V+ +v-)]^ ~ 1.839^, for fc = 3, 

where = (l9 ± 3 \/33 )^^^ . 

From the above calculations, note that the number of non-jammed sectors put together, i.e. the sum of all those 
with lengths £ G [1, L — 1], increases as fast as the number of sectors with C — L. Since these latter can not evolve any 
further, each constitutes a separate sector having only one configuration (as opposed to non-jammed strings which, 

from the ASEP analogy, bear (^) state configurations). 

Growth of sectors with Nq = 0.- Following this line of reasoning, it is straightforward to determine also the number 
of IS conserving the internal energy throughout. That is the situation referred to after Eq. ([7|), where no Aq characters 
appear in the IS. For this sector we now define Gc (1) and Gc (0) as the number of invariant C- strings having no 
consecutive O's, and whose first bit is 1 and respectively. Thus, the counting of strings constrained by Nq — 
requires the evaluation of Gc (1)- Clearly, by construction these numbers involve the relations 

Gc{Q) = Gc-i(l), 

Gcil) = Gc{Ai) + ... + GciAk-i), (17) 

where, as before, GciAi) refers to the number of irreducible sectors, now subject to Nq = 0, having Ai as their first 
character (1 < i < fc — 1). Also, it can be readily verified that these latter numbers are involved recursively in the 
same form as their _F's counterparts in Eq. p^. When using those relations in Eq. p7)) . the following linear recurrence 
immediately emerges 

fc 

G£(1) = ^G£_,(1) , Ok. (18) 

1=2 

The characteris tic p olynomials Pk{x) associated to the generic term of this sequence now distinguish the parity of 
fc-mers, namely [24| 

X]i=o x^ — x^, for fc odd , 
Pk{x) = { (19) 
X]i=o a;^* — a;*^"-^, for fc even. 
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which along with the boundary terms Gi(l),... ,Gk{l) determine the specific form of Gc V£ > fc. As expected, 
the roots of the above polynomials only yield exponential growth for k > 2 (evidently for dimers G/;(l) = 1, in 
correspondence with the sole [Ai]^^^ configuration). In the limit £ — > cx) these sectors grow progressively faster as 
k increases but in all cases slower than the respective F's of the unrestricted sectors [Eq. ([Ti)) ]. as they should. In 
particular, for trimers it turns out that 

G£(l) cx ~ 1.325^, fc = 3, (20) 



where w 




III. NUMERICAL RESULTS 



Armed with the ASEP correspondence discussed before, we have conducted extensive simulations of SS currents and 
pair correlations in several subspaces for both dimers and trimers. In all cases, we evolved indpcndcnt configurations 
for each of the studied sectors. The corresponding initial conditions were prepared by random deposition of {L — C)/k 
ASEP monomers, that is A^. particles (fc = 2, 3), on a ring of Nq + A^i + • • • + Nf^ sites. Here, we distinguished 
between fc — 1 type of ASEP vacancies, so we tagged them in the same order as that appearing in the particular 
IS's considered, either periodic or not. Subsequently, each ASEP particle was duplicated (triplicated for the case of 
trimers), by adding another particle (two particles) over one (two) extra adjacent location(s) specially created for 
that purpose, i.e. 1 — >• 11 (1 111). In turn, the tagged vacancies were replaced accordingly by j consecutive I's 
followed by a 0, that is — >■ 1 • • • 1 0, if they referred to Aj characters, while keeping all ^o's as O's. This defines an 
efficient algorithm to produce generic configurations in the chosen IS sector within the original DRLG representation 
of L sites. Owing to PBC there is a small hindrance however, as eventually cyclic shifts in one site (one or two, 
for trimers) might be necessary to maintain invariable all sublattice densities in the generated DRLG configurations. 
These latter were then updated with the stochastic rules summarized in Table HI using chains of L = 1.2 x 10* sites 
evolving typically up to 10^ simulation steps. Each of these ones involved L update attempts at random locations, 
after which the time scale was increased by one unit, i.e. t — >■ t + 1, irrespective of these attempts being successful. 

The above algorithm enabled us to average measurements over nearly 5 x 10^ histories originated from independent 
sector configurations, thus reducing significantly the scatter of our data. We considered three typical periodic situations 
which were afterwards compared with the null string [ analogous to the monomer DLG, as already mentioned in Eqs. ^ 
[7])]- These are specified in Table HIl along with their sublattice saturation currents and densities, in turn arising from 
Eq. (jlOp and simple stochiometric considerations. We also examined random strings generated from the ASEP version 
of these periodic sectors by swapping randomly through the lattice the order of their irreducible characters, thus 
keeping all relative concentrations. 



A. Currents 



The measurement of SS sublattice currents in these sectors involved the monitoring of transient regimes which 
for most temperatures and fields decayed typically in ^ 10'^ steps. Then, we averaged all currents along two 
further decades during which no significant fluctuations were observed. As usual, the sublattice currents J7i — 
( — N~ )^ j^^^ can be defined operationally using the total number of forwards and backwards N~ par- 
ticle jumps within sublattice A^, and averaging over all event realizations during an interval {t, t + At). However to 
avoid any dependence on that latter lapse, particularly inconvenient to monitor early non-stationary stages, instead 
we measured instantaneous correlators of the form 



Jt{Ek,l3J,t) — ( (nj+i . . . Hj+k-i) {Rj n.j fij+k - Rj n-j+k rij) )^ 



1 1 

- ± - tanh 

2 2 



Eh 



+ 2/3J (rij+fc+i - rij^i) 



(21) 



where fij = 1 — rij are just vacancy occupations numbers in sublattice A^. Here, ( )t denotes an ensemble average over 
these correlators at time i, whereas right (left) hoppings i?"*" (i?~) are defined so as to take into account the rates 
referred to in Table H] In Fig.[T] we show the resulting SS currents normalized to the saturation values of Table [Hi 
after taking T/ J = ±1, ±0.1, and C/L = 1/2 for several driving fields. In addition. Fig. [5] displays other plane of the 
current, in this case holding = 10 for a variety of temperatures, using both attractive and repulsive couplings. 
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FIG. 1: (Color online) Normalized SS currents for reconstructing trimers and dimers (insets) using LjL = 1/2. In listed order, 
the sectors of Table lU here correspond to squares, triangles, and circles respectively. In (a) T j J = 1, and (b) Tj J = 0.1, all 
data collapse onto the null string current (for clarity not shown), upon rescaling the drifts Ek with a's such that |1 — a| ^ 0.05. 
The data are fitted by the mean field currents of Refs. [i,|23l (solid lines), using the associated DLG densities referred to in 
Eq. ((8]). The mean field description extends as well to (c) T/J = —1, and (d) T/J = —0.1, where however all currents are 
sector dependent and DLG densities are chosen only as fitting parameters. 
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J: IS sector Jsat / tanh {^^) Density 

2 [AiA^f'' '"-'i^'r'^ P=l-t 



2 lAA,f^' (apzii^ P=l-i§ 

2 lAAlf/' 



(l-pi)(l-2pi) 

Pi -2 ( pi = 1 



(1~P2) (1-2P2) 
P2-2 

2P3 (1-P3) 
3-P3 



'2L 



2P2 (1-P2) i p2 = 1 - § 

3-P2 

A A 1^/** (p-i)(2p-i) n - 1 _L 




f4 A 431^/8 (p-l)(8p-3) „_1 5£ 

[^2^1^oJ 7p-12 P — J- 8L 



TABLE II: Limiting sublattice currents and densities for dimers and trimers {k = 2, 3) in periodic strings of length C considered 
in the simulations below. These are formed by repeating, e.g. [Ai^o] = [(10)(10)(0)] C/5 times, etc. In some sectors these 
quantities are common to all sublattices. 



It turns out that normalized currents of nonequivalent sublattices are indistinguishable within our error margins, a 
non obvious feature (except for random strings, as their sublattice occupations approach each other as £ — > oo). 

(i) J > 0.- More importantly, under attractive interactions the currents of all studied sectors, both periodic and 
random, can be made to collapse into that of the null string by rescaling slightly the driving fields. This is displayed 
in Figs.[T^ and[T]D (random sectors not shown, for legibility). Furthermore, the data collapse extends also to the J — T 
plane provided the attractive couplings are taken slimly rescaled, as illustrated in Fig. [2^. It is noteworthy that in 
both J — Ek and J ~ T planes the null sector data follow very closely the mean field currents of the usual DLG 
pol | . These arise essentially from a kinetic version of the cluster variation method applied to dynamics proceeding via 
exchange processes As expected from the arguments given in Sec. II A, here the fitting of the null string current 
is attained upon choosing — (l — x) / [-'^ + r('^ ~ ^ ^^'^ monomer density for the DLG system [see Eq. (|8]) ]. 
These numerical observations naturally lead us to put forward the universality hypothesis 

for normalized sublattice currents in ferromagnetic DRLG . Here, / — f{PJ) and g = g{Ek) are sector dependent 
scaling factors (probably close to 1), whereas J^]^^ is taken as in Eq. pH)) . Preliminary runs using other string lengths 
indicate similar results, thus adding more weight to this conjecture. 

(ii) J < 0.- On the other hand, for repulsive interactions all normalized currents come out to be sector dependent, 
as is shown in Figs.[lJ;, [TJl, and [2b, an aspect becoming more pronounced as temperature is lowered. However, this 
dependence appears to involve only the relative concentration of irreducible characters rather than their particular 
distributions, because in all situations the normalized currents of random sectors follow closely those of their periodic 
counterparts (alike the attractive case). Although these currents can also be fitted using the mean field approach to 
DLG UllOl, the corresponding monomer densities can no longer be understood in terms of the ASEP analogy given 
before (except for null sectors and/or high drive regimes, as expected). Here, we merely use DLG densities as fitting 
parameters which actually turn out to be rather sensitive in tuning the values of the current plateaux exhibited in 
Figs.[Ui, and [lb. 

Further to universality issues, or the lack thereof for J < 0, these results also suggest that DRLG currents inherit 
some of the salient features of DLG ones according to the type of interactions [1, HSl (cf- nevertheless pair correlations 
in non null sectors). For J > there is a continuous changeover from a rather insulating state at low temperatures 
and fields to a conducting phase as both T and £ are increased. For J < however, already at low temperatures 
some sectors can be found in conducting states, while exhibiting comparatively larger conductivities at small fields. 
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0.5 1 0.5 1 

y T/j -T/j 

FIG. 2: (Color online) Normalized SS currents versus temperature for fc = 3, and k — 2 (insets) while keeping £k/T = 10 for 
(a) J > 0, and (b) J < 0. Key to symbols is the same as in Fig.[T] In (a) the data collapse is obtained by the rescale J — >■ J/7, 
choosing 7's such that |1 — 7I ^ 0.02. These currents follow closely the null string data (not displayed to avoid overcrowding). 
By contrast, in (b) all currents are sector dependent. In both cases the mean field fittings (Refs. [J, [20|, solid lines) arise from 
the considerations given for Fig.[T] 

As for the appearance of current plateaux in Figs.[T}i, andt^b, notice that these are already present at the mean field 
level of the standard DLG. 



B. Pair correlations 

Turning to mesoscopic scales, in the following we focus on the SS instantaneous density-density correlation functions 
expressed through the subtracted or cummulant form 

'^(^^ = ;^ XI ( ^ ) ~ J , (23) 

j 

with p\. being the density of sublattice A^. Also, to gain some further insight into the average organization of 
stationary regimes, we consider the static structure factor or Fourier transform of C(r) 

5(g) = i ^''''C{r). (24) 

-L/2 < r < L/2 

which in our case is a real function of the wavelength A = In j q. In Fig.[3]we first show this latter function for the case 
of dimers and trimers in the null sector, taking £/L — 1/3 and 1/4 in turn. At finite temperatures and far from the 
saturation or ASEP regime \E\ ^ | J|, the main maxima might be regarded as remnants of the periodicity and long 
range order of the Ising ground states. Clearly, in the attractive and repulsive situations these are respectively of the 
form yl^/('''+^)^^/('^+^' and [A^. ^o]^^'*^^^' (along with translations), so as T — !• it is natural to expect sharp peaks 
at g = and in each case. However, in general notice that when C/L ^ the ground state of the null sector 
is highly degenerate for J < [having the form A^^Aq A^'^Aq ■ ■ ■ with mj > constrained as J2j "^j — ~ 
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thus setting a residual entropy which grows linearly with the system size. In fact, as temperature is lowered the 
structure factor remains essentially broad (see data of T/J = —1/2 for ZjL — 1/2), and correlation lengths are of 
the order of the lattice spacing. A similar scenario arises in the monomer DLG, where degeneracies for J < and 
p 7^ 1/2 also preclude long range order at any temperature 



1.5 



S(q) 

1 



0.5 





71/2 71 

q 

FIG. 3: (Color online) Structure factors of sector \Ai^'~ (null string) for dimers and trimers (inset), using CIL = 1/3 and 
1/4 respectively. Solid, dashed, and dotted lines stand in turn for (Ek,T / J) — (2, — 1),(2, 1), and (8, ±1). The two first 
cases exhibit maxima very near to g = 27r/3 (dimers), n/2 (trimers), and (both), thus resembling the average equilibrium 
periodicity of low temperature regimes. At high fields no periodicity takes over, regardless of the sign of J. The dashed-dotted 
lines denote the situation {Ek,T/J) = (2, —1/2) for £/L = 1/2, where the ground state is highly degenerate. 

(i) J < - Yet, a rather different situation shows up for repulsive couplings in other periodic sectors also bearing 
high degeneracies. This is observed in Figs.|4^ and |41d where, as an example, we illustrate respectively the behavior 
of strings [A^Ao]^/^ (fc = 2), and [A2AiAl]^/^ {k = 3). It can be readily verified that for J < the number of 
ground states of the former grows exponentially with the lattice size so long as C/L ^ 5/7, whereas that of the latter 
grows in the same manner provided C/L ^ 8/17 (otherwise these states would be plain periodic sequences of the 
form [A^yl2^o]^^^ for dimers, and [A2A1 {A3 Aq)'^]-'"^^'^ for trimers). However, despite that exponential degeneracy 
the structure factors of both sectors can single out wavelengths exhibiting peaks that rapidly narrow and heighten as 
temperature is lowered. Also in lowering T, correlation lengths turn out to grow monotonically (e.g. at T/J = —1/2, 
^ ~ 500 in sector [A2 Ai Aq]^^^), while being already relatively large at T/J = — 1 as compared with those of the 
case J > (see below) . This is suggestive of an order-by-disorder scenario [21[ in which thermal fluctuations are able 
to lift part of the degeneracy by selecting a subset of states with largest entropy in the ground state manifold. Here, 
note that the role of frustration in Refs. [21| is being played by the string conservation laws and the spatial extent of 
characters. 

(ii) J > 0.- By contrast, under attractive interactions no ordering seems to emerge for these sectors at low 
temperature regimes. Structure factors now remain basically fiat (i.e. bounded) and correlation lengths do not 
increase as T is lowered (Figs. 2^ andSb)- Thus, thermal fiuctuations now appear as being unable to suppress the 
residual entropy (also growing linearly in the thermodynamic limit of both sectors V £/i 7^ 0'*', 1~) and no order-by- 
disorder seems likely to occur, just as in the situation of the null string and standard DLG under repulsive couplings. 
Similar results for J > were observed in the other sectors of Table [Hi instead exhibiting effects of order by thermal 
fiuctuations as long as the couplings are taken repulsive. 

(Hi) \£\ ^ \ J\~ In approaching the ASEP or high field limit, pair correlations become progressively independent 
of the coupling values because all configurations tend to be equally likely. Already for {Ek,T/J) = (8,±1), Figs.|31 
13^, andHJj indicate that all results are numerically indistinguishable. However as the field is lowered from the ASEP 
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FIG. 4: (Color online) Stationary structure factors for (a) dimers in sector [yljAp]'''''^, (b) trimers in [Aj Ag]''"''*, (c) trimers 
in energy conserving sectors [A^A J''"''* (solid line, circles), and [^jAJ^'''' (dashed line, squares), using LjL = 1/2 throughout. 
Spatial correlation functions (insets) are highly oscillating, but only local maxima evidencing a decay oc e~^^^ are shown. 
Although (a) and (b) bear both highly degenerate ground states irrespective of the coupling sign, for J < the peaks marked 
by arrows rapidly grow as temperature is lowered. As in Fig.O here solid, dashed, and dotted lines (in turn, circles, squares, 
and triangles for insets) denote respectively the cases [Ek,T/J) = (2, — 1), (2, 1), and (8, ±1). In the high field regime the 
results of (a) and (b) become independent of T/J, whereas in (c) this is always the case whether or not the system is driven. 
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regime, for J < correlations can enhance significantly their ranges, in turn becoming arbitrarily long if temperatures 
are taken low enough. As displayed in Figs.|3^ and|3)D, this is to be contrasted with the opposite trend of both <^ and 
S{q) under attractive couplings, where the drive favors a slight increase of pair correlations. 

(iv) Nq = 0.- Finally, the ASEP regime is also related to the energy conserving sectors (A^o — 0) referred to after 
Eq. (O, for which the presence of NN interactions is inconsequential. Despite the fact that every state is equally 
weighted, owing to the fc-mer size here structure factors can still single out characteristic wave numbers (sector 
dependent) and exhibit large correlation lengths. This is illustrated in Fig.[3J;. In the ASEP limit of the null sector 
these issues were recently analyzed in Ref. ,18j where closed expressions for general /c's were given for static correlations. 
Such rich behavior contrasts to that of totally uncorrelated monomers and dimers in sectors [Ai]^^'^ (the only ones 
with A'o = 0, which can be viewed as hopping monomers within one independent sublattice). 

IV. CONCLUDING REMARKS 

To summarize, we have discussed stationary aspects of driven lattice gases in which the role of biased monomers 
is played by extended and reconstructing fc-mers under a Kawasaki dynamics [ Eq. 1^ and Table U] . Exploiting the 
correspondence between these processes and those involving the particle species defined in Eq. ^ we readily identified 
the many sector decomposition of the original problem, ultimately encompassed in the invariant ordering of these 
characters along the so called irreducible strings [l^ [l^ . In the high field regime the dynamics of these new particles 
was thought of as an asymmetric exclusion process defined on a smaller lattice (Sec. II A), thereby enabling us to 
evaluate saturation currents [Eq. (|10[ ) ] for generic strings or sectors of motion. In turn, the proliferation of these 
latter was shown to grow exponentially with the length of these strings (Sec. II B). 

At finite temperatures and drives, we studied numerically the case of both dimers and trimers in typical sectors 
whose initial configurations were prepared by random sequential adsorption of 'monomers' in the equivalent ASEP 
states. These were then raised up into DRLG configurations, always keeping the distribution of irreducible characters 
or tagged ASEP 'vacancies' in each of the studied sectors (Table |ll|. The emerging stationary currents clearly 
discern between universal (Figs.[T^,[T}D,[5^), and sector dependent behavior (Figs.[Tt,[TJi, according to the particle 
couplings being attractive or repulsive. In the former case, the universality hypothesis put forward in Eq. ([^ suggests 
in turn an effective medium relation between generic sectors and null string currents via a slight rescale of drifts and 
interactions. 

When it comes to mesoscopic levels of description (Sec. IIIB), also distinctive features appeared at low temperature 
regimes. In spite of the residual entropy in most of the studied sectors (stemming from their highly degenerate ground 
states), under repulsive couplings there is a substantial increase of both correlation lengths and structure factors at 
characteristic wavenumbers as temperature is lowered (Figs.S^ andjUa). We interpret these modes as being selected 
by thermal fluctuations from the ground state manifold thus giving rise, we suggest, to an order-by-disorder scenario 
[2lj. For attractive interactions however, this latter can not be inferred from our simulations since, in line with the 
ground states degeneracies, both ^ and S{q) do not grow any further; a situation which resembles that of the standard 
DLG {p ^ 1/2), and null strings (£/i ^ '^pi) under repulsive couplings (Fig. [3]). Finally, both large drive regimes 
and non-interacting sectors are governed by a simple product measure, but due to the fc-mer size their correlation 
functions may show nontrivial oscillations and large correlation lengths (Figs.|4^, |4}d, and|4j;) [l^. 

It is natural to ask whether the above numerical findings could be approached theoretically. At the microscopic 
level of the master equation j23 | , the formal analogies between this latter and the Schrodinger equation describing the 
evolution of associated quantum spin chains have proven useful in the analysis of several nonequilibrium processes 
[(3, 26]. In fact, for vanishing drives and interactions these reconstructing systems have been studied in terms of spin-i 
Heisenberg ferromagnets [l3|, but for £, /3J ^ the evolution operators are neither familiar nor simple to analyze. 
On the other hand, already at the mean field level it is not clear how to proceed with an exponential number of 
conservation laws such as the IS's discussed throughout. 

Other issues not covered here that would be worth pursuing concern the phase ordering dynamics of periodic 
sectors under repulsive couplings where, as we have pointed out, stationary correlation lengths can get very large 
at low temperature regimes. It would be intersting to determine whether the dynamic exponents characterizing the 
large time growth of f actually depend on the subspaces where the evolution takes place. There is also the question 
about tagged particle diffusion either with or without driving fields. For /3J = and £ ^ it is known that in 
the non-reconstructing case the root-mean-square displacement of a tagged particle around its mean position grows 
asymptotically in time as ~ t^^^ (as usual), but if the bias is zero it grows anomalously slow as ~ t^^^ [13]. In the 
reconstructing situation, where such caging effect might depend on the particular distribution of fragments of the 
sector considered, these issues remain quite open. The question also extends to the coarsening regime — /3 J ^ 1 for 
which other displacement laws might emerge depending on whether or not the motion is driven [28| . 
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